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Abstract 



with real 



An elementary set of non-Hermitian N by N matrices (A) 7^ H^ N > (A) 
spectra is considered, assuming that each of these matrices represents a selfadjoint 
quantum Hamiltonian in an ad hoc Hilbert space of states n iphysical) . The problem 
of an explicit specification of all of these spaces (i.e., in essence, of all of the eligible 
ad hoc inner products and metric operators 0) is addressed. The problem is shown 
exactly solvable and, for every size N — 2, 4, . . . and parameter A G (—1, 1) in matrix 
H^ N \\), the complete iV— parametric set of metrics 0^ a2 ... a (A) is recurrently 
defined by closed formula. 



1 Introduction 



1.1 The concept of the metric 6 in Hilbert space 7-^0% SJCa O 

One-dimensional quantum systems described, in units h = 2m = 1, by the ordinary 
differential Schrodinger equation 

d 2 

HxP(x) =Ei){x), H = -— + V{x) (1) 

serve as a universal testing ground for the ideas, methods and techniques of quantum 
mechanics. One works with the standard representation L 2 (M) of the Hilbert space 
of states where the bound states are normalized in usual manner and where the 
Hamiltonian itself is self-adjoint, 

/ i/;*(x)ip(x)dx = l, H = H ] . (2) 

JR 

In parallel, the scattering solutions ip(x) of eq. ([1]) offer the simplest illustration of 
the delocalized waves which must remain compatible with the unitarity of the time 
evolution, etc. (cf., e.g., ref. [T] for numerous illustrations). 

The transparency of such an elementary implementation of quantum theory can 
prove deceptive. People often forget that the requirement of the Hermiticity of H in 
L 2 (M) can be replaced by an alternative, equally acceptable requirement H = of 
the Hermiticity of the same operator in another Hilbert space f{ip h y sical ) ^ L 2 (]R) 
where a different definition of the inner product (ip, tp') e is employed, 

(V>,<0')e= / ^*{x)Q{x,x')^\x')dxdx' , H = H x = Q^H^Q. (3) 

Although the space L 2 (M) remains unchanged as a vector space, the definition of the 
"correct" linear functionals (i.e., the definition of the mapping T of the "ket- vectors" 
upon "bra- vectors" ) is less elementary in ^(physical) ^ 

Tip{x)= I xp*(x')0(x',x)dx' . (4) 

JR 

Misunderstandings may emerge whenever the metric is nontrivial, 9 = 0^/. 
Then, our Hamiltonian H = H$ appears non-Hermitian in the conventional Hilbert 
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space f{( Dirac ) — L 2 (]R). The paradox has an elementary resolution since the space 
n (Dirac) with itg stanc i ar d 

inner product is not the space of states of the system 
in question (cf. Appendix A for a brief recollection of a few concrete illustrative 
examples with this "cryptohermiticity" [2] property). 

1.2 The problem of the ambiguity of the metric G = 0(H) 

The deeper study of the similarity ([3]) between an operator H and its adjoint W in 
a preselected space Ji^ Dtrac ) dates back to the early sixties 0j. In physics, the first 
use of such a feature, i.e., of the so called quasi- Hermiticity constraint 

OH = H^O (5) 

imposed upon a sufficiently nontrivial and realistic Hamiltonian H ^ W emerged 
much later [3]. Still, up to now the subject remains full of open questions. One 
of the most challenging ones concerns the ambiguity of the metric O assigned, via 
eq. (j5J), to a given Hamiltonian operator H. Indeed, eq. (JSJ) itself defines "too many" 
alternative physical metric operators = 0(H). An explicit constructive illustration 
of such an ambiguity of the assignment H — > Q given in section 5 of ref. [5] employs 
a free Hamiltonian H = H which is complemented by a two-parametric family of 
non-trivial metric operators 

Q (M ostaf azadeh) ^ K)=e ~F _ <p^hK) (@) 

where V denotes parity and where F and K are arbitrary real numbers. 

The first discussion of the general problem of the ambiguity of 0(H) has been 
published by Scholtz et al [3]. They emphasized that besides the Hamiltonian H 
itself, any other operator O = Oj of an observable quantity in ^}{(p h v sical ) must obey 
the same Hermiticity relation as H = Oq. In an opposite direction, any eligible 
physical metric operator Q must remain compatible with the corresponding set of 
the quasi-Hermiticity relations 

eOi = oje, j = (o),i,2,.... (7) 
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These requirements reduce the ambiguity in at every index j. In this sense, the 
choice of the physical metric = 0(H) can be made, in principle, unique. 

We intend to return to the problem of the ambiguity of the metric in what 
follows. One of our reasons is that the universal strategy represented by requirement 
(jTJ) is in fact rarely successful in practice. The solution of the complete set of the 
linear operator relations ([Z]) appears to be hardly feasible. Typically, just j = 
is considered and a particular solution 0(H) of eq. (jSJ) is sought for. Appendix 
B reviews a few alternative proposals of making the Hamiltonian-dependent metric 
operator 0(H) unique in such a case. 

1.3 Unique metrics in certain matrix models of scattering 

A schematic comparison of a few alternative techniques of the removal of the am- 
biguity of metrics 0(H) ^ I has been performed in our two brief comments [5J. 
We restricted our attention to the mere two-dimensional Hilbert spaces r \-^P h y stcal ) _ 
Using an elementary set of two-by-two matrices H = H^ we compared the merits 
and shortcomings of various versions of non-equivalent 0s. Due to the simplicity of 
the space we were able to base our analysis on an explicit construction of all the 
solutions 0^ of eq. flS}. 

In some sense we shall just extend such a study to certain less trivial Hamiltonian 
matrices H = H^ in what follows. The practical feasibility of such a project relies 
on suitable simplifications. It is obvious that for a general matrix H( N ' one could 
hardly consider its size in the range N > 4 [7j. 

One of the most natural, anharmonic-oscillator-inspired choices of the simplified 
tridiagonal matrices H^ was proposed in ref. [SJ. Even these models with the 
number of variable matrix elements limited to N/2 appeared to be only tractable 
numerically [9j. In ref. [ID] this observation led us to the most drastic reduction 
of the allowed number k of the variable matrix elements in H^ N '(Xi, . . . , A&) to the 
smallest integers k = 1, 2 and 3. 

It was a pleasant surprise to discover that the latter choice proved extraordinarily 
successful. Without difficulties we were able to consider all the matrix sizes N = 2K 
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including even the limiting case of N = oo. The "exceptional" choice of N = oo 
made us ready to study and solve certain difficult conceptual problems in scattering 
theory (cf. the text of ref. [TU] for more details). One of the reasons was that already 
the one-parametric matrix model of dynamics 



H = H(X) 



2 -1 






-1 2 


-1 




-1 


2 


-1 - A 




-1 + A 


2 -1 






-1 2 -1 






-1 2 '•. 



(8) 



proved compatible with the diagonal matrix solution of eq. fl5J), 



e[F(A)] 



1-A 




1 - A 




1-A 






1 + A 




1 + A 




1 + A 



(9) 



The existence of such a local metric was already considered improbable in the ded- 
icated phenomenological literature [H] and we felt encouraged us to generalize the 
explicit formula (Q to Hamiltonians with more parameters (cf. [10]) and to some 
matrices with different structure (cf. [T2]). 

Here, due to the lack of space, we shall skip all the similar enhancements of so- 
phistication. Rather, we shall return to the bound-state problems where N < oo in 
what follows. Before doing so we should add a remark that there exists an amaz- 
ingly close and direct connection between the two apparently independent sample 



Hamiltonians given by eqs. ([I]) and (jHJ). In Appendices C ad D we shall explain this 
relationship in more detail. 



1.4 Matrix models and bound states 



In our present return to model (jHJ) and to the bound states we shall consider all the 
truncated, finite-dimensional matrix descendants H^ N ' of eq. (jHJ) with truncations 
N = IK. For this sequence of one-parametric toy Hamiltonians 



H^(X) 



H (4 \X) 








H^(X) 
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1 + A 

-1 
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-1 + A 




-1 
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-1 + A 





-1 - A 
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-1 — A 

2 -1 

-1 2 



(10) 



'11} 






1 - A 
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-1 








-1 
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-1 










'121 



etc the main result of our present paper will be the explicit construction of the 
respective complete sets of all the metrics Q^ N '(X) in closed form. 

Such a project is nontrivial since the bound-state wave functions remain localized 
so that there is no point in demanding the asymptotic locality constraint which made 
the scattering metric virtually unique in [10] . One can arrive at a really satisfactory 
physical interpretation of the bound-state system only via an exhaustive knowledge 
of all the metrics Q(H) allowed by eq. (jSJ). 

In an introductory part of our present paper we shall set A = in H^ N \X). In 
section [2] this simplification will help us to explain our method of construction of all 
the admissible metrics 0^(0) = Oq^- In essence, we shall combine the brute-force 
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symbolic-manipulation constructions performed at the first few N = 2,4,... with 
the subsequent extrapolation of the resulting closed formulae towards all the even 
integers N = 2K. 

In the second half of our paper (cf. sections [3] and H]) we shall return to the 
nontrivial, asymmetric Hamiltonian matrices H = H^ N \X) with A 7^ 0. Firstly, 
in section [3j we shall solve eq. ([5]) for the first three models ( fTUj) . ( fTTj) and (fl2l) . 
In section H] we shall then extrapolate the resulting triplet of metrics to all the 
superscripts N = 2K. The closed formula for all of the solutions Q( N >(\) of eq. (J5]) 
will be obtained as our main result. 

In section we shall summarize our message while Appendices A - D will com- 
plement it by a few additional remarks and technical notes. 

2 The description of the method: A = 



2.1 Trivial starting point: All the metrics at N = 2. 

In the light of ref. |6J the simplest possible two-dimensional Hamiltonian 



2 -1 
-1 2 



is easily assigned the real pair of energies E± = 2 ± 1 as well as the general real-matrix 
ansatz for the metric 

b r 

reflecting its necessary Hermiticity. In an encouraging start of our systematic study 
we may insert both these matrices in eq. (jSJ) and arrive at the single constraint / = /'. 

All the resulting two-parametric metrics 0o possess eigenvalues 9± expressible in 
closed form, 6± = f ±b. It is trivial to conclude that our O is positive (and can be 
called a metric) iff / = /' > and f 2 > b 2 . 
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2.2 All the metrics 6 at N = 4. 



In the first nontrivial step of our analysis let us consider the four-dimensional Hamil- 
tonian H^(X) at A = 0, 



H, 



(4) 



ijW(O) 



2-100 
- 1 2-1 
0-1 2-1 
0-12 

and let us try to deduce the generic form of all of the related matrices 0o directly from 
the set of iV 2 = 16 equations (jSJ) for the N 2 = 16 unknown (though, presumably, 
real) matrix elements of Go- These equations are not all linearly independent. No 
surprise - the general solution contains N = 4 real parameters [TBI [H] . 



2.2.1 Construction 

As a typical task for Mathematica or Maple we solved our set by the brute force 
methods of linear algebra and we obtained its complete four-parametric solution 

«2 OL\ + «3 «2 + «4 Ois 
«3 «2 + «4 Oil + as Ctg 
«4 OJ3 «2 Oil 

exhibiting linear dependence on all of its four parameters, 



9> v 



(4) 



(13) 







(i) 



ai Mi + a 2 M 2 + a 3 M 3 + a 4 M 4 . 



(14) 



While is just the four-dimensional unit matrix, the remaining three expansion 
matrices represent its elementary sparse-matrix generalizations, 



M. 



("1) 



10 
10 10 
10 1 
10 



M, 



("1) 



M 



(4) 



10 
10 1 
10 10 
10 

This result of the computation indicates the possibility of the existence of a certain 
friendly extrapolation pattern towards the metrics Oq^ at any higher N. 



1 
10 
10 
10 
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2.2.2 Positivity 

After we specify the Hamiltonian but before we select any particular solution = 
0(H) = ©"I" of eq. (j^D we have to guarantee that our candidate for the metric is 
invertible and positive definite. Only then, this operator can consistently specify the 
corresponding physical Hilbert space 

^physical) of stateg of 

our quantum system [I]. 

At the larger dimensions N the proof of the positivity may be difficult. In the 
four-dimensional matrix example ffTSl) it degenerates to the mere four elementary 
inequalities 

-2 a 4 + 2 a x + a 3 - a 2 ± Vh (—a 2 + "3) > , 

2 a 4 + 2 ot\ + «3 + «g ± \fh («2 + CK3) > . 

They must be satisfied as a guarantee of the positivity of all the four eigenvalues 8k 
of the metric O^ . 

2.3 Extrapolation to N > 4 
2.3.1 An ansatz for 6^° 

It is natural to expect that formula (I14p is just the first special case of the general 
expansion 

Oi N) = J2 aj M^(0). (15) 
i=i 

Let us activate the experience collected at the smallest N and assume that all of 
the matrices M = Mj (0) are solely composed of the matrix elements or 1. In 
the j—th matrix the location of all of the non- vanishing elements may tentatively be 
selected as follows, 

(itff ) (0) = 1 iff i-k = m, N + l-i-k = n, (16) 

m = j- l,j-Z,...,l-j, n = N - j, N - j - 2, . . . , j - N . 

Such an educated guess generalizes the above N = 2 and N = 4 results to all even 
dimensions. Its validity has carefully been verified at several higher even integers 
N = 2K. One should note that the mere insertion of the ansatz followed by the 
check of the result is quick. 

9 



2.3.2 Verification: N = 6 etc. 



B (6) 



Formulae (fT5l) and ffTBl) determine all the extrapolated 2K— parametric matrices 
Q( 2K ). For illustrative purposes let us pick up N = 2K = 6. This choice gives the 
formula 

a.} a 2 ct$ ol^ a 5 a@ 

a 2 a,i + ol 3 a 2 + a 4 a 3 + a 5 a 4 + a 6 a 5 

a 3 a 2 + OL Jt a} + a 3 + a 5 a 2 + a 4 + a 6 a 3 + a 5 a 4 
a 4 a 3 + a 5 a 2 + a 4 + a 6 + a 3 + a 5 a 2 + a 4 a 3 
a 5 a 4 + a 6 a 3 + a 5 a 2 + a 4 + a 3 a 2 

olq 0L5 a 4 a 3 a 2 oti 

We may easily verify the validity of the pertaining set of linear equations §5§ by the 
simple-minded and straightforward insertion again. 

A different category of verification of the internal consistency of our general 
N— parametric N < 00 result f|T5|) results from a direct study of the continuums 
limit iV — > 00. A few comments on this interesting are collected in Appendix D. 



3 Metrics 0^ for H = H^ N \\) at A > and N <6 
3.1 Model with N = 2 

The key features of bound states in our discrete short-range interaction models 
become already well illustrated via their most elementary special case (fTUj) . Firstly, 
this is the simplest model which shares peculiarity of the spectra which remain real in 
the N— independent interval of couplings A G (—1,1). For the whole sequence of our 
Hamiltonians we shall parametrize A = cos ip G (—1, 1), therefore, with cp G (0, tt/2). 

At iV = 2 we can easily evaluate not only the closed formula for the energies, 
E = E±^ = 2 ± sin<y2, but also the norm = 2 sin 2 ip of the related eigenstates ip±. 
In addition, the cryptohermitian model H^ 2 '(cos(p) nicely illustrates the difference 
between its right eigenvectors and their left-eigenstate partners. In the spirit and no- 
tation of ref. [H] these respective column vectors \ip±) and their row-vector partners 
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are different, 



IV>±) 



1 + COS if \ I 1 — COS 09 

, IV±» ~ 

=F sin Lp I \ =F sin <p 

but a biorthogonal basis can be formed of them. Thus, H ^ (cos cp) is a self-adjoint 
matrix in an ad hoc, Hamiltonian-dependent Hilbert space of states / }{(P h y s ' lcal ) . 

In the light of refs. [6] a key merit of the N = 2 example can be seen in the 
existence of the explicit spectral definition of the metric, 



e = HM>*+«iM +IV>-»*- 



(17) 



In this representation the guarantee of the necessary positivity of the metric reads 
t± > 0. After the insertion of the eigenvectors we arrive at another explicit formula 
for the metric, 
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(1 — cosip) 2 {t + + t_) (1 — cosy?) sin</?(— 1 + + tJ 
[1 — cos <p) sin (/?(— 1+ + t_) sin 2 <p(t + + t_) 

Its inspection reveals that the metric may be re-written as a superposition 
6 (2) (A) = an m{ 2 \\) + a 2 M 2 (2) (A) 



(19) 



with the two A— dependent matrix coefficients, 



M X (2) (A) 



1 - A 
1 + A 



, M 2 (2) (A) 



1 

1 



(20) 



Such a re-parametrization leaves the positivity criterion entirely transparent, 



e*i > 0, a 2 (l - A 2 ) > c^, N = 2 



(21) 



so that we may choose any a 2 from the interval (— a\ sin a?, a.\ sinyj). The transition 
to the Hermitian limit A — > appears facilitated in the new parametrization. 



3.2 Shorthand notation 

The continuity of the expansion matrices in the free-motion limit as noticed above 
remains true at all N > 2. Thus, we may visualize metrics 0^(A) as expanded in 
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terms of the generalized A— dependent sparse matrix coefficients obtained as a certain 
A— deformation of their A = predecessors defined by closed formula (fT6|) . With this 
perspective in mind let us now define the following infinite sequence of polynomials, 

P = l, Pl ±) = l±\, P 2 = l-A 2 , P 3 (±) = (1 ± A) (1 - A 2 ) , 

P 4 = (l-A 2 ) 2 , P 5 {±) = (1±A)(1-A 2 ) 2 , P 6 = (l-A 2 ) 3 , .... (22) 
In terms of these polynomials the doublet of our sparse expansion matrices M^lW 

(2) 

can be characterized by the "incidence" or "indexing" matrices S{ 2 with certain 
integer (or empty) entries. In general they will carry all the information about the 
position and about the degree of polynomial matrix elements P n = P n (A) of the 
respective matrices Mj N \\). At N = 2 they are defined simply by the following 
assignment, 

1 

1 

1 



In what follows we shall demonstrate, step by step, that the polynomials P^+^A) 
with the minus-sign superscript will always sit in the left upper triangle (i.e., above 
the second diagonal) of the respective expansion coefficient Mj (A) and vice versa. 
Thus, at any N, the expansion-coefficient functions Mj N \\) will be unambiguously 
determined by the mere indexing matrices Sj N K Of course, up to now we only 
demonstrated that these observations are valid at N = 2. 

3.3 Model with N = 4 

Hamiltonian H^ N \X) of eq. (ITTj) is the simplest model which is purely kinetic near 
its "distant" lattice points and which is dynamically nontrivial just in the vicinity 
of the origin. The A— dependent coupling merely connects the two points xn/2 and 
xn/2+i in the middle of the lattice. The four eigenvalues of matrix H^(X) remain 



m x (2) (a) 



PI 



(-) 



pi 



(+) 



(2) 



Mf (A) 






Po 




Po 
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Figure 1: Spectrum of H^(X). 
real in the same interval of couplings A G (—1, 1) as above (cf. Figured]), 



£±,± = 2 ± - V6 - 2 A 2 ± 2 V5 - 6 A 2 + A 4 . 

Symbolic manipulations on the computer enable us to find all the corresponding 
matrices of the metric 0^(A), 

(1 - A) a 2 (1 - A) a 3 a 4 

a 2 (1 — A) oci (1 - A) + oc 3 (1 - A) o^+a 2 (l-A 2 ) a 5 

a 3 a^+a a (l-A 2 ) (1 + A) + a 3 (1 + A) a 2 (1 + A) 

«4 a 5 a 2 (1 + A) ax (1 + A) 

They may again be interpreted as the sums 

6 (4) (A) = on Mi + a 2 M 2 +a 3 M 3 + a A M A (23) 

where 



Mi 



1-A 

1-A 

1+A 

1+A 



Mo 



1-A 

1-A 1-A 2 

1-A 2 1 + A 

1+A 



13 



10 

1 - A 1 

1 1+A 
1 



Ma 



1 
10 
10 
10 



(24) 



In the shorthand notation of our previous paragraph the following four incidence 
matrices Sj will carry again all the necessary information about the respective four 
matrix polynomial functions of A. In computations, these incidence matrices S*- , 



i.e., 



1 

1 2 
2 1 
1 





1 
1 








will be used for the encoding and/or efficient reconstruction of the respective expan- 
sion matrices Mj (A). 

It is worth noticing that even the simplest metric with a 2 = «3 = 0J4 which 
is proportional to the first coefficient Mi and which remains diagonal (i.e., in the 
language of coordinates on the lattice, "local" ) ceases to be proportional to the unit 
matrix so that our model resides in a nontrivial Hilbert space where 9 / /. 



3.4 Model with N = 6 

Although all the six eigenvalues of the matrix if( 6 )(A) may be expressed in closed 
form in principle, we shall only graphically confirm that all of them remain real in the 
same interval as above, with A G (—1,1) (cf. Figure [2]). Inside this interval the metric 
0( 6 )(A) exists and its general form is obtainable from eq. ([SD by its straightforward 
computer-assisted solution. The resulting matrices 0^ 6 ^(A) are displayed here in the 
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Figure 2: Spectrum of # (6) (A). 

following split presentation, 

a! (1 — A) a 2 (1 — A) a 3 (1 — A) 

0^(1 -A) ctj (1 - A) + a 3 (1 - A) a 2 (1 - A) + a 4 (1 - A) 

a 5 (l-A) a 2 (l-A) + e^(l-A) a, (1 - A) + a 5 (1 - A) (1 - A 2 ) + a 5 (1 - A) ... 

a 4 a 3 (l-\ 2 ) + a 5 a 2 (1 - A 2 ) + a 4 (1 - A 2 ) + a 6 

a 5 a 4 + a 6 a 3 (1 — A 2 ) + a 5 

a 6 a 5 a 4 

a 4 a 5 a 6 

a 3 (1 — A 2 ) + a 5 a 4 +a 6 a 5 

a 2 (1 - A 2 ) + a 4 (1 - A 2 ) + a 6 a 3 (1 - A 2 ) + a 5 a 4 

... a! (l + A) + a 5 (l + A)(l-A 2 ) + « 5 (l + A) a 2 (1 + A) + a 4 (1 + A) a 3 (1 + A) 

a 2 (l + \) + a 4 (1 + A) a, (1 + A) + a 3 (1 + A) <xg (1 + A) 

a 5 (1 + A) a 2 (1 + A) a, (1 + A) 

The use of the shorthand symbols becomes indispensable for the sufficiently 
efficient and compact encoding of the N = 6 expansion formula 



N 



e w (A) = E«i< , W 

i=i 



(25) 
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At j — 1, 2, . . . , 6 the six respective compact incidence or indexing matrices Sf^ are 







1 1 
1 2 










2 1 
1 1 










1 2 

1 3 2 

2 3 1 
2 1 

1 







These matrices form an inseparable part of any practical computer-assisted applica- 
tion of the formalism at N > 4. At N = 6, moreover, their explicit form also offers 
an interesting insight in their N— dependence forming a sufficiently inspiring starting 
point of our extrapolation programme. 



4 Extrapolation 

Starting from N = 8, the explicit form of the indexing matrices Sj 1 *'* becomes also 
rather large for being printed. Still, their computer-assisted use remains as easy and 
straightforward as at N = 10 etc. Thus, in the main part of our homework we simply 
formulated and tested the alternative extrapolation hypotheses. Now, it remains for 
us just to summarize the results. 



16 



4.1 Recurrences for the incidence matrices, with the central 
ones exempted. 

It is easy to reconstruct the A— dependent matrices Mj 2K \\) from the knowledge 
of the respective shorthand symbols Sj 2K \ What remains for us to do is to define 
all the set of the shorthand symbols Sj representing all our expansion matrices. 
In such a context, the trial and error method enabled us to collect a sufficiently 
extensive set of these symbols. In the next step, their inspection revealed that at 
any given N = 2K the first K — 1 matrices Sf K ^ with j = 1,2, . . . , K — 1 would be 
easily constructed from their predecessors Sf K . 

Purely empirically, the latter recurrent construction has been found in an en- 
largement of the dimension followed by a symmetric attachment of the two j — plets 
of units "1" in the empty parts of the left upper corner and of the right lower corner. 

In an entirely similar manner, the last K matrices S^/f-i-i-j w ^h j = 1,2, ... ,K 
become also formed in the similar manner. Explicitly, their K predecessors S^k-x-j 
must be modified by attaching j zeros "0" in the right upper corner and in the left 
lower corner. 

In both of the "leftmost-subsequence" and "rightmost-subsequence" scenarios, 
the results displayed in section [3] offer a sufficiently instructive illustration of such 
a recipe. An explicit algebraic reformulation of such a doublet of two-dimensional 
recurrences would be also as straightforward and compact as their above, purely 
verbal description. An algorithmic version of these recurrences was, after all, needed 
also during our practical computer-assisted evaluation of the matrix elements of 
& N >{X) at all the N which we considered. Of course, all these descriptions of our 
two-dimensional recurrences are equivalent and amenable to the rigorous proof by 
mathematical induction. Due to the lack of space, the routine details of algebra as 
well as of its rigorous proofs are skipped here and left to interested readers. 
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4.2 Recurrences for the central incidence matrices S 



(2K) 
K 



At the remaining subscript j — K the construction of the most complicated missing 
member S^ K ^ of the family must be discussed separately. In both the "leftmost- 
subsequence" and "rightmost-subsequence" scenarios its (2K — 2) by (2K — 2) pre- 
decessor proves, rather unexpectedly different from the naively expected matrix 

Q (2K-2) 

This means that the sequence of the "middle" or "central" matrices should 
be treated as exceptional. Fortunately, they remain created by a straightforward 
recurrent recipe. Its idea relies on the use of certain specific predecessor matri- 
ces CS 2K ~ 2 \ With a freedom in their specification let us decide to proceed in the 
"rightmost-subsequence" manner. This means we shall enlarge the dimension of C 
() whatever it is) and we shall fill K "neighboring" units "1" in the left upper corner 
and in the right lower corner. 

We are now ready to define the specific predecessors C( 2K ~ 2 \ They appear to be 
constructed from the old "middle" matrices via a specific two-step recipe. 

Firstly we replace each "old" numerical element in S^^ 2 ^ by its successor, i.e., we 
replace "old 0" by "1", "old 1" by "2", etc. In the second step we form a left- 
right reflection of the resulting matrix and arrive at the final form of the necessary 
predecessor CS 2K - 2 ^ as a result. Thus, at N = 4 we have the sequence 



s{ 2) = 



£(2) = 



Similarly, the recurrent construction of the matrix at N = 6 will result from 
adding six units "1" to the auxiliary predecessor matrix in the formula 



2 

2 3 
3 2 
2 



£(4) 



2 

3 2 
2 3 
2 



etc. We may conclude that the central matrices S^ K ^ at the respective K — 1, 2, 3, 4 
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(etc) form the sequence 

1 

1 2 

1 3 2 
13 4 2 

2 4 3 1 

2 3 1 
2 1 

1 

(etc). The general pattern of their recurrent construction is obvious. 
4.3 Verifications 

The full-fledged formulae for the eight-parametric 0^ 8 ^(A) already cease to be easily 
printable but their characterization using the incidence matrices remains fully trans- 
parent and compact. All of the individual expansion matrices entering the general 
series (1251) for metrics exhibit the same simultaneous change of the sign of A after 
the reflection with respect to their second diagonal. This is well visible in our last 
illustrative equation 

1-A 

1 - A 1-A 2 
1-A (1-A) (1-A 2 ) 

1-A (1-A) (1-A 2 ) (1-A 2 ) 2 

1-A 2 (1-A 2 ) 2 

1-A 2 (1 + A) (1-A 2 ) ... 

1-A 2 

1 + A 

where a part of the real and symmetric matrix M^ K \\) at K = 4 is displayed. 



1 

1 2 
2 1 
1 



1 

1 2 

1 3 2 

2 3 1 
2 1 

1 
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5 Summary 



In quantum theory an operator H represents an observable provided only that it is 
self-adjoint in a Hilbert space equipped with a metric 9. For a given H, equation 
H^Q = QH specifies a complete menu of all the eligible = 0(H) needed to 
determine the inner product. We illustrated the feasibility of the construction of 
all of these Gs for an infinite sequence of certain one-parametric N by N matrices 
ifW(A) with A G (-1,1). 

A recurrent method of the construction of all the admissible metric matrices 
9^(A) has been proposed and tested at A = 0. For A G (—1, 1), the straightforward 
construction of the individual N— parametric @^(A) = 0^ Q2 ajv (A) was based 
on the computer-assisted symbolic manipulations at the smallest iV < N( minima iy 
Next we extrapolated these formulae for the metric to all the subsequent larger 
N > N( minima iy We carefully verified the validity of our extrapolations via tests 
performed at a number of sample N > N( minima iy 

We revealed that the recurrences which are needed for the reconstruction of all 
the set of all the iV— parametric matrices O^(A) = 0^p a2 a N W fr° m the set of 
their N — 1 predecessors ( ^ 7V ~ 1 ^(A) = < 3> ( ^^ t l]... :CtN _ 1 (^) degenerate to the recurrences 
needed for the reconstruction of the iV— plets of matrix coefficients Mj N \\). The 
ultimate simplification of our recurrent recipe has been achieved when all the matrices 
Mj N \\) proved easily constructed from the knowledge of the related elementary 
indexing matrices Sj N ^ with integer or empty entries. 

During our study of our matrix toy-model Hamiltonians we persuaded ourselves 
that the underlying mathematics is friendly and that the recurrences generate the 
matrix indices Sj N ^ in a really transparent form. Although this fact should par- 
tially be attributed to the mere one-parametric choice of our class of Hamiltonians 
H( N \\), we firmly believe that our method which mixed the low— N evaluations 
with subsequent all— iV extrapolations might keep its efficiency for a number of more 
realistic and, in particular, more-parametric sparse-matrix toy-model Hamiltonians. 
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Figure captions 

Figure 1. Spectrum of H^(X). 
Figure 2. Spectrum of H^(X). 
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Appendix A: Toy-model potentials in eq. (QD) 

Among the oldest explicit examples of the phenomenological Hamiltonian operators 
H exhibiting the cryptohermiticity property Q one may cite, e.g., the imaginary 
cubic anharmonic oscillators studied by Caliceti et al [15], the "wrong-sign" quartic 
anharmonic oscillators described by Buslaev and Grecchi and others [16] or a family 
of the more general, "non- Hermit ian" but stable bound-state models as proposed 
by Bender and Boettcher [17J. In all of these one- dimensional schematic models 
of bound states the reconstruction of any suitable metric from a given oscillator 
Hamiltonian H proves rather difficult. Pars pro toto we may recall the perturbation- 
series construction of 0(H) for the imaginary cubic H [5], with many further relevant 
references cited therein. 

A perceivably better picture of the properties of the physics-determining metrics 
0(H) has been obtained for certain exactly solvable potentials of bound states (cf. 
the semi- numerical construction of 0(H) for a square- well model H ^ W [18] as an 
example) . 

The situation appeared perceivably worsened after transition to the scattering 
regime where even the combination of a sophisticated perturbation construction with 
the choice of the really most elementary exactly solvable interaction potentials in 
eq. (JT]) did leave many conceptual questions unanswered (cf., e.g., refs. [TT| fT9]). 

Appendix B: Conventional choices of = 0(H) in 
quantum theory 

In the majority of applications of quantum theory using non-Dirac Hilbert spaces 
^(physical) construction of the metric remains almost prohibitively difficult even 
when one restricts attention to the single observable, i.e., to j — in eq. ([7]). One 
may feel forced to work just with a drastically simplified and/or highly schematic 
class, say, of point-interaction Hamiltonians [191 [20]. The construction of 0(H) may 
remain feasible only when one employs an approximate (e.g., perturbation [5J [TB"] ) 
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method. Still, even in many papers which define a quantized system in an unusual 
Hilbert space f{(p h y sical ) which is equipped with a non-Dirac metric / J their 
authors usually avoid the overcomplicated requirement ([7]) and select just one of 
particular solutions Q(H) of eq. §5§ instead. 

As one of the the simplest illustrations of such a strategy we might recall even 
the most common model (CQ) with a most common real potential. Evidently, just the 
trivial solution 6 {Dirac) = / of eq. ([5]) is being assigned to H = H^. In particular, 
out of all of the above-cited eligible metrics Q^ Mosta f azadeh ^p^ ^ j us ^ m ost trivial 
particular solution with F = K = is being used in connection with the free-motion 
version of eq. ([I]). 

In VT— symmetric quantum mechanics admitting B ^ Q{D%rac) thorough 
recent review written by Carl Bender [21]), the problem of the ambiguity of the choice 
of the metric O = Q(H) has been circumvented as well. Although this less traditional 
formalism admits various nonstandard, apparently non-Hermitian models (including 
even field models with real spectra |22j etc), the current choice of the space / }{ ( J> h y stcal ) 
is equally restrictive, preferring special metrics Q( Bender ) — QV where C = C(H) 
represents a unique "charge" while V is the usual parity. 

The most natural generalization of the VT— symmetric theories with G ^ Q(Dirac) 
has been described by Mostafazadeh [13J. He re-attracted the attention of the in- 
ternational scientific community to the abstract quantization rule (j7]) of ref. [1] and 
to the related ambiguity of the reconstruction of the correct Hilbert space f{(p h y sical ) 
form a given Hamiltonian. He worked out some illustrative examples (cf. [5|) and, 
together with Batal [18], he emphasized the possible physical relevance of metrics 

Q -J- Q (-Bender) 



25 



Appendix C: Discretized Runge-Kutta version of 



eq. (ffl) 

One of the key simplifications of some of the technical aspects of solving differen- 
tial Schrodinger eq. (CQ) is commonly sought in its replacement by its Runge-Kutta 
difference-equation approximation 



*(*™>-2*M + *(*^> + v(Xk) iixt) = Ei , (Xk) 



(26) 



(cf., e.g., ref. [23] for more details). In place of the real line of coordinates x G M 
the equidistant lattice of points x^ may be conveniently defined by the formula 
Xk = — 1 + kh, k G Z in terms of a suitable (i.e., usually, sufficiently small) real 
constant h > 0. For the purposes of the description of scattering the lattice remains 
infinite so that the kinetic energy operator —d 2 /dx 2 may be visualized as the following 
tridiagonal matrix 



Hn 



2 
-1 



-1 

2 -1 
-I 2 



(27) 



— I 

2 -1 

— I 2 



In the bound-state context with the Dirichlet boundary conditions ip(— 1) = ^(1) = 
one usually considers just the finite set of the lattice points, k = 1, 2, . . . , N. With 
Xn+i = +1 we, in effect, fix an elementary length h = 2/(N + 1) which would vanish 
in the continuum limit N — ► oo. 

In the latter scenario the explicit specification of the non-Dirac Hilbert space 
^(physical) mus ^ slightly modified, 



(2* 
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On this background the numerical use of the approximation ([I]) — > (126]) finds numer- 
ous applications in cryptohermitian quantum mechanics. In papers [21], for example, 
the exact, analytic solvability of the differential eq. ([T]) with a certain class of suf- 
ficiently simple complex potentials V(x) 7^ V*(x) has been shown paralleled by the 
exact solvability of the discrete partner eq. (1261) • 

Whenever the potential V(x) remains sufficiently smooth, the role of the dis- 
cretization errors may be expected negligible, in the domain of the sufficiently large 
iV ^> 1 at least. Difficulties may arise for point interactions. They proved popular 
[25] and found applications in relativistic equations [26] and in manybody systems 
[27] . This, as we already mentioned, motivated our interest in in the new field of 
cryptounitary scattering [T2] and, in particular, in its simplest discrete model ([S]). 

Appendix D: Models H^ N \\) at large N > 1 

At a sufficiently large N our one-parametric Hamiltonian H^ N '(X) represents in fact 
one of the discrete versions of eq. (CEJ) with a certain point interaction V(x) localized 
in the origin. For a deeper understanding of such a correspondence let us asccept 
that xq = — 1 and x 2 k+i = 1 and let us abbreviate 2 — h 2 E = cose as usual [T2| |2"3]. 
We may then treat e G (0, 7r) as a new energy variable and we may visualize the wave 
functions ip(x) with i ^ as satisfying the free-motion eq. (I26I) . complemented by 
the boundary conditions ip(xo) = i/j(x2k+i) = 0. This picture of the N ^ 1 system 
must be completed by the doublet of the A— dependent relations 

(1 + A) i/)(x K+1 ) - 2 cos e ip(x K ) + ^(xjc-i) = , (29) 

^{x K+2 ) - 2cose^(xx+i) + (1 - \)ip(x K ) = (30) 

so that we may assume the emergence of a discontinuity of ip{x) in the origin. 

At the sufficiently large N ~ 1/h ^> 1 the wave functions near x = remain 
well represented by their respective one-sided Taylor series so that eqs. (|29l and (13"0|) 
may be interpreted simply as a matching condition. We return to the original energy 
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variable h 2 E = 2 — 2 cose = F and insert the truncated expansions 

Mk-i) = Mo) - + o{h 2 ) , Mk) = mo) - \ hMo) + o(h 2 ) , 

Mk+i) = Mo) + ~ hMo) + o(h 2 ) , Mk+2) = Mo) + \hMo) + o(h 2 ) 

in eqs. fl29|) and fl30l) . A straightforward algebra leads to the following elementary 
condition 



h 
2 





Mo) 



(31) 



\ MO) j \f-i i-\ J \ mo) 

which matches the wave functions and their derivatives in the origin. In the domain 
of sufficiently small h > this relation is equivalent to the original constraints fT29|) 
and fl30|) . We may conclude that in the continuum limit iV — > oo our sequence 
of the matrix Hamiltonians H^ N \\) can be reinterpreted as a series of dynamical 
models which converge to a specific differential eq. ([T]) which is split in two halves. 
Indeed, the h — > limit of eq. (1511) produces the elementary opaque-wall constraint 

Mo) = Mo) = o. 

At all the nonvanishing small "elementary lengths" h > 0, our rigorous defini- 
tion fl3Tl) leaves the point-interaction potential term V(a;) translucent and manifestly 
energy- dependent. Its definition by the mixed boundary conditions is nonlinear in the 
coupling A. Various special cases of this N ^ 1 bound-state model may be studied 
noticing, for example, that the energy-dependence disappears in the low-excitation 
regime where the quantity F = h 2 E remains negligible. 

Another interesting special case is encountered when the interaction is completely 
switched off, A — ► 0. Then we may recollect our iV— parametric free-motion formula 
(|15p for the metric 0o and we may compare it with its two-parametric differential- 
operator counterpart ([6]). This comparison confirms that a linear combination of the 
matrices M\ (0) and Mjy (0) survives the limiting transition iV —>■ oo while, in 
contrast, all the other terms in (|15|) disappear due to the build-in requirement of the 
absence of an elementary length in the Mostafazadeh's theory of ref. [5]. 
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